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Dark matter density profiles based upon ACDM cosmology motivate an ansatz velocity distribu- 
tion function with fewer high velocity particles than the Maxwell-Boltzmann distribution or proposed 
variants. The high velocity tail of the distribution is determined by the outer slope of the dark mat- 
ter halo - the large radius behavior of the Galactic dark matter density. N-body simulations of 
Galactic halos reproduce the high velocity behavior of this ansatz. Predictions for direct detection 
rates are dramatically affected for models where the threshold scattering velocity is within 30% of 
the escape velocity. 
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I. INTRODUCTION 

Dark matter is the dominant form of matter in the 
Universe and measuring its properties is one of the most 
important outstanding problems in astrophysics and par- 
ticle physics. Directly detecting Galactic dark matter 
through terrestrial experiments is probably the least am- 
biguous method for determining its identity. Accurately 
predicting scattering rates for direct detection experi- 
ments requires a model for the velocity distribution of 
Galactic dark matter. This article proposes a class of ve- 
locity distribution functions that are motivated by cos- 
mological models for the Galactic dark matter and stud- 
ies how this ansatz affects predictions for direct detection 
experiments. 

Dark matter halos form through a continuous process 
of smooth accretion and merging of smaller mass halos 
in the ACDM model of cosmology. Clues about the for- 
mation history of these halos can be found in their phase 
space distributions, motivating a careful study of dark 
matter densities and velocities. Modern numerical sim- 
ulations, which probe a large range of halo mass scales, 
have shown that the merging process leads to a near uni- 
versal double-power-law density profile that is approx- 
imately described by the Navarro-Frenk- White (NFW) 
model [1 . Though the physical origin for NFW-like pro- 
files is unknown, they provide a valuable link between 
the small scale properties of dark matter halos and the 
cosmological model on large scales. 

The velocity distribution is dynamically related to the 
density profile for equilibrium dark matter halos - In 
direct dark matter detection studies, the 'Standard Halo 
Model' (SHM) 4] is the canonical velocity distribution 
for interpreting results from experiments and constrain- 
ing dark matter models [5] . Isothermal and isotropic ap- 
proximations to the dark matter distribution functions 
such as the SHM are useful because they simplify rate 
calculations. However, recent high resolution cosmolog- 
ical simulations provide strong evidence that the SHM 
does not appropriately capture the behavior of the dark 
matter in the Galactic halo [BHS]. By mapping out the 
dark matter phase space distribution, these N-body sim- 
ulations consistently point to deviations from a standard 



Maxwellian velocity distribution function, most promi- 
nently at the tail of the distribution 10J. This new ev- 
idence motivates a reevaluation of the SHM and careful 
consideration of velocity distributions that better capture 
the features of the Galactic halo. 

There have been recent attempts to derive the dark 
matter velocity distribution from first principles or 
from observations of rotation curves [H]. However, the 
resulting models over-predict the number of high veloc- 
ity particles compared to the highest resolution dark 
matter-only numerical simulations with ~ 10 3 Mq parti- 
cles [HI [7] , as well as simulations that include baryons [5] . 
It has been claimed that Tsallis distributions, derived 
from non-extensive statistical mechanics |13j . provide a 
good fit for the baryon plus dark matter simulation in [5] . 
However, because the resolution for this simulation is 
~ 7 x 10 5 Mq, it is limited at probing the high veloc- 
ity tail relative to the dark matter only simulations. An 
alternate approach is to fit the numerical results with 
an arbitrary function jTJ [9] ; while these models are use- 
ful for making experimental predictions, there is no clear 
motivation for the fit parameters. 

This paper obtains a velocity distribution function that 
describes the double power-law density models. While 
these velocity distributions have been obtained numeri- 
cally for isotropic models and also for certain classes of 
spherical anisotropic models [TlUl6j . no analytic solu- 
tion has been found. The isotropic velocity distribution 
function corresponding to double power-law density pro- 
files PQ is described by 
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where k is the power-law index, v csc is the escape veloc- 
ity, and Vq is the dispersion. This distribution function is 
significantly different from Maxwell-Boltzmann distribu- 
tions and its variants in the high velocity tail. It describes 
only the smooth dark matter component; modifications 
to the distribution function from discrete subhalos and 
streams can be important phenomenologically but need 
to be added separately [I?Ml9| . Additionally, deviations 
from a smooth distribution arise from gravitational scat- 
tering from the Sun and the Jovian planets |20j . 
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Direct dark matter detection experiments look for low 
energy nuclear recoil events caused by dark matter scat- 
tering off of target nuclei [3T] . The nuclear recoil thresh- 
olds for dark matter detection experiments are in the 
range of tens of keV, which implies that the recoil ener- 
gies in many models of dark matter are well below thresh- 
old. Only particles with the largest velocities have suf- 
ficiently energetic nuclear recoils to be visible at direct 
detection experiments. This fact makes it particularly 
important to accurately describe the high velocity tail 
of the dark matter velocity distribution in the halo sur- 
rounding the Milky Way. Relative to a Maxwell distribu- 
tion, predicts a smaller fraction of dark matter near 
the escape velocity, which affects the relative rates be- 
tween direct detection experiments, especially those with 
different target nuclei or energy thresholds. 

This paper is organized as follows. Sec. |ll] motivates 
the ansatz in ([l ) and derives a general form for the power- 
law index. Sec. Ill compares to numerical calculations 
of the full velocity distribution function for Milky Way- 
like spatial profiles, including a gravitational potential 
model of the disk. The ansatz in ([I]) is also compared 
to the results of N-body simulations. Elementary direct 
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detection phenomenology is explored in Sec. IV and con- 
clusions are presented in Sec. |Vj 



II. EQUILIBRIUM DARK MATTER 
DISTRIBUTIONS 

The Galactic halo forms through a process of hier- 
archical merging, smooth accretion, and violent relax- 
ation [U |52 [53] ■ The central regions of halos, where 
relaxation processes have ended, are in quasi-static equi- 
librium |24j . In this case the dominant component of 
Galactic dark matter is described by a steady-state dis- 
tribution function that evolves slowly with time. In the 
outer regions of the halo, at low binding energies, the 
phase space distribution will depend more dramatically 
on both the history and environment of the halo. For 
example, a halo with a fairly quiescent formation history 
will have a different phase space distribution at low bind- 
ing energies compared to one with a recent major merger 
[5J. Sec. Ill argues that the behavior of the distribu- 



tion function of the Milky Way at high velocities and low 
binding energies resembles an equilibrated system on av- 
erage. There can be large excursions from equilibrium 
at the high velocity tail due to recent accretion; how- 
ever, N-body simulations suggest that these are spatially 
localized streams. 

According to Jeans theorem, the six-dimensional phase 
space distribution function f(x,v) for a spherically sym- 
metric and isotropic system can be written in terms of 
the energy, which is an integral of motion: 



£ oc 2iP(x) - v 2 



(2) 



£ oc v csc (r) 



(3) 



Due to the assumption of equilibrium, /(£) for £ < 
must vanish since this corresponds to particles with ve- 
locities greater than the escape velocity of the system. 

If the dark matter of the Milky Way is nearly equi- 
librated, then the Jeans theorem implies that the dark 
matter density determines the velocity distribution func- 
tion. Dark matter halos in cosmological simulations are 
well fit by spherically averaged double power-law density 
distributions of the form, 



p(r) 
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(4) 



where p s is the scale density, r s is the scale radius, a 
is the slope of the halo density near the core, and 7 is 
the slope at large radii. For example the NFW profile 
has (a, 7) = (1,3) pQ, while the Hernquist model has 
(a, 7) = (1,4) [15]. The Jaffc model has (a, 7) = (2,4), 
but the a = 2 inner slope is steeper than indicated by 
numerical simulations. 

Velocity distribution functions that behave as 
Maxwell-Boltzmann distribution functions at low veloc- 
ities arise from dark matter halos with inner slopes of 
a = 2. This article is primarily concerned with the ef- 
fects of the dark matter velocity distribution on direct 
detection rates. The relative velocity between dark mat- 
ter and target nuclei for v <C v is set by the solar veloc- 
ity; therefore, the low velocity behavior of the Galactic 
dark matter does not significantly affect direct detection 
rates. Though finding a velocity distribution function 
that describes a < 2 is important, these modifications 
to the dark matter velocity do not qualitatively change 
the direct detection rates [5B] . The goal of this section is 
to derive a phenomenological velocity distribution func- 
tion that can arise from these double power-law density 
models to predict direct dark matter detection rates. 

A given density distribution is related to the velocity 
distribution through the gravitational potential of the 
halo, ip(r): 



V 2 V = -4ttGp = -AttGM / f(v)d 3 v 



(5) 



where G is the gravitational constant and M is the halo 
mass. For an isotropic density distribution, there is a one- 
to-one correspondence between the density and ergodic 
distribution function, given by the Eddington formula: 
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The local escape velocity is defined as the velocity where 



The expression for f(ip) can be written in terms of the 
binding energy £ by replacing ijj — > £ . Equation ^ 
is a powerful simplification that allows one to solve for 
the distribution function of an arbitrary spherical density 
model. In many cases, however, the solution for f(v) is 
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analytically intractable and must be obtained by solving 
|6]) numerically. This is the case for most combinations 
of (a, 7) in Q; as a result, there is no general closed- 
form analytic solution for the velocity distribution func- 
tion that corresponds to double power-law density mod- 
els. Fortunately, it is still possible to assume an ansatz 
for f(v) that reproduces the numerical solution of the 
Eddington formula for double power-law densities. The 
possibilities for such an ansatz are limited if one assumes 
that the dark matter halo is in equilibrium. 

The high velocity tail of the local dark matter distribu- 
tion function arising from equilibrated, double power-law 
models is in conflict with the Standard Halo Model and 
its variants. Consider the following generalization of the 
Standard Halo Model as an ansatz 



f(£) ex (e £ / £ ° - l) k Q(£), 



(7) 



where k is the power-law index of the distribution and 
describes the behavior near the escape velocity. The ben- 
efit of this ansatz is that it satisfies the Jeans theorem for 
an equilibrated system and goes continuously to zero at 
the escape velocity. 

The power-law index is defined as 



k = lim k(£ ), 



with 



k(£) 



£ df{£) 
f(£) d£ ■ 



(8) 



(9) 



For double-power density profiles of the form in Q , the 
power-law index can be evaluated analytically. In this 
case, the second term of the Eddington formula is negli- 
gible. Expanding the density, potential, and the distri- 
bution function in ^ around small £ gives a power-law 
index 



k = 7 — 



(10) 



for 7 > 3 [57]. As 7 — > 3, k(£) does not approach a simple 
power-law and the approximations that lead to ( 10 1 break 
down 1141. Sec. Ill fits the numerical solutions to the 



parameters for the ansatz velocity distribution function 
in ([!]). As 7 — > 3, the best fits for k tend to be k ~ 



2.0, slightly larger than (10). Earlier studies from galaxy 
formation via violent relaxation motivated k — 1.5 128 . 



It is remarkable that k takes such a general and simple 
form. The index is determined almost exclusively by the 
outer slope of the density distribution; all terms that de- 
pend on a vanish in the low energy limit. The outer slope 
controls the behavior of k because the dark matter parti- 
cles with the highest velocities are those with the small- 
est binding energies. These particles will be in highly 
energetic orbits about the halo, and will be concentrated 
at large radii, far from the core. Density distributions 
with larger outer slopes have fewer particles orbiting at 
large radii, which means that the low-energy component 



of f(£) is suppressed. A large value of k precisely cap- 
tures this behavior. 

The distribution function for binding energies can be 
rewritten in terms of velocities using the relation in ^ 
to get (JTJ . The velocity distribution ansatz in (JlJ is well- 
described by a Gaussian peaked near vq for v <§; v csc - As 
v — > v csc , the distribution function approaches 



f(v) -> (v csc - v) k . 



(11) 



Cosmological N-body simulations indicate 7^3 — 5 [1 , 
|2T)] . which means that the velocity distribution falls off 
near the escape velocity to the power k — [1.5,3.5]. A 
more detailed comparison with cosmological simulations 
follows in Sec. IIII1 

The Standard Halo Model and the King model provide 
a useful comparison to the distribution in |7| and they 
are defined, respectively, as 



/shm(£) = N(£ a )e £ / £ °0(£) 
/ Ki ng(£) - N{£v)(e £ / £ ° - 1)Q(S). 



(12) 



These distributions are frequently used for direct dark 
matter detection predictions because they make the cal- 
culations tractable. In addition, they satisfy Jeans theo- 
rem under the assumption of isotropy and spherical sym- 
metry. However, they do not correspond to NFW-like 
density models, especially near the high-velocity tails of 
the distributions. In particular, the SHM behaves near 
the tail as k — > 0, and the King model has k = 1. Conse- 
quently, these velocity distribution functions over-predict 
the number of particles in the tail of the distribution. 

The Tsallis distribution is another model for the veloc- 
ity distribution that has been recently discussed in the 
literature, and is defined as 



9/(1-9) 



/TsaUis(v) OC 1 - (1 - q)- 



(13) 



The Tsallis distribution predicts that the escape velocity 
is given by Wg SC = Vq/(1 — q) and k = q/(l — q), where 
the distribution in (JlJ sets these three parameters inde- 
pendently. Hansen et al. [3D] show there is a correlation 
between the parameter q and the local density slope, im- 
plying that q varies with radius if the density slope is 
not the same at all radii. A disadvantage of this model, 
however, is that does not satisfy the Jeans theorem for 
spherical and isotropic systems if the circular velocity, v , 
is held constant. This violation of Jeans theorem is also 
true for the generalized Maxwellian distribution, which 
has been used to model the radial and tangential compo- 
nents of the velocity distribution [SJ. 

The formalism and models in this section only apply 
to the spherically-averaged velocity distribution, and do 
not capture any physics pertaining to steams, subhalos, 
or any other structure in the phase space distribution. 
Though streams are unlikely to affect the overall velocity 
distribution because their densities are less than ~ 0.1% 
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of the smooth dark matter density [BJ, they are seen 
to occupy the high energy tail of the velocity distribu- 
tion mus]. 



III. NUMERICAL EVALUATION 

The previous section derived how the high velocity 
tail of the dark matter phase space is populated for 
cosmologically-motivated spatial densities. The ansatz in 
([I]) interpolates between a Maxwellian distribution at low 
velocities to one matching the high velocity behavior in 
(JlOfc. This section compares the analytic approximation 



in ( |10[ ) with the velocity distributions that are calculated 
from the full numerical solution of the Eddington formula 
in ^ . Additionally, it shows that the ansatz distribution 
function accurately describes the velocity distributions of 
simulated dark matter halos both with and without bary- 
onic physics. 



A. Distribution Function Model for the Milky Way 

A full gravitational potential model for the Milky Way 
is necessary to evaluate the Eddington formula, which 
requires including a model for the baryons. The potential 
model is comprised of the sum of three components: the 
dark matter, the disk, and the bulge, so that Vtot = 
"0haio + V'disk + V'buigc- The halo mass within a given 
radius is obtained from Q, while the velocity of the local 
standard of rest (LSR) and the escape velocity are given 
by 



dipt, 



J \ST 



dr 



(14) 



respectively. 

The value of f/'haio is uniquely determined for a given 
pair, (a, 7), and for specific values of the scale radius, r s , 
and scale density, p s . The velocity distribution ansatz 
in ([T]) is compatible with any value of the outer slope 
with 7 > 3, while the inner slope is given by a = 2. 
However, the behavior of the potential for r<r s dom- 
inantly determines the low velocity behavior of Galactic 
dark matter and does not significantly affect predictions 
for direct detection rates. 

The full model for '(/'total needs to include additional 
potentials for the bulge and disk components. Here, 
a spherically-symmetric potential is used for the bulge 
component of the form 



V'buigc = -GM huLse /r, 



(15) 



where the mass of the bulge is Mb u i gc = 10 10 M [3j. The 
disk is modeled by a spherical potential of the form 



disk 



-GM dWt (l-e- p / r -)/r. 



(16) 



The potential formula in (16) contributes ~ 115 km s 
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FIG. 1: Region of parameter space in local density (defined 
at 8.5 kpc) and scale radius, r s , consistent with (i) halo mass 
internal to 60 kpc in the range (4±0.7)xl0 11 M Q , (ii) escape 
velocity in the range 450-600 km s _1 , and (iii) circular velocity 
in the range 200-280 km s" 1 , for models with a = 1 and 7 = 3 
(red), 4 (yellow), and 5 (green). 



stars with M^k = 5 x 10 M and a scale length of 
= 5 kpc. The spherical approximation in ( 16 ) is 



the monopole term for a standard axisymmetric double 
exponential disk potential [3J. In comparison, for the as- 
sumed disk mass, the axisymmetric double exponential 
disk potential gives v\ sr ~ 130 km s _1 at the solar radius, 
and is 11% larger than vi SI derived from the spherical ap- 
proximation. The disk potential sets the circular velocity 
of the Milky Way that arises from non-dark matter, i.e. 
visible matter. 

The measured values of i>i sr , v csc , and the integrated 
halo mass provide constraints on the parameters of the 
velocity distribution function. The circular velocity is in 
the range [33] 



Uta = [200,280] km s~ x . 



(17) 



The 90% confidence level on the escape velocity is in the 
range 



[498,608] km s" 1 , 



(18) 



and is determined from local high velocity stars by the 
RAVE survey [33J . This constraint on the escape velocity 
is only attainable for a prior on the spectral index for the 
velocity distribution of halo stars, 



with 



oc (v csc -v) k % 



K = [2.7,4.7]. 



(19) 



(20) 



Simulations indicate that the power-law index for halo 
stars is steeper than dark matter. The values of £;* in 
(20) correspond to k — [0.5,2.5] for dark matter [33| . 



to the circular velocity for a disk consisting of gas and The lower values of k correspond to unphysical halos, 
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FIG. 2: Distribution functions for models with a = 1 and 
7 = 3 (red), 4 (yellow), and 5 (green) obtained by numerically 
solving the Eddington equation with a disk potential as in 
(16 1. The dashed lines have slopes of 7 — 3/2. The bottom 
panel shows the error of using the ansatz velocity distribution 
(solid) or the King model (dotted) instead of the complete 
numerical solution. 



with fc* near the upper end of the range in ( 20 ) . Since 



the escape velocity is correlated with the spectral index 
of the stars and the dark matter, if values of k and fc* on 
the lower end of the allowed region are used, lower values 
of v esc must be used, i.e. v csc ~ [450, 550] km s _1 . 

Finally, there is a bound on the integrated dark matter 
mass of the Milky Way halo within 60 kpc, determined 
by distant halo stars to be [31] 



M{< 60 kpc) = (4.0 ±0.7) x 10 U M, 



(21) 



In the ACDM model, an extrapolation of the mass profile 
to larger radii implies a total Milky Way halo mass of 
~ 10 12 M [331|35]. The above constraints are examined 
in the context of the power-law parameterization in Q; 
more generic halo profiles are beyond the scope of this 
analysis, see [35] for further discussion. 

Fig. [l] shows the regime of scale radius-local dark mat- 
ter density parameter space that is consistent with the 
observational bounds above, where the local dark matter 
density is simply defined as p(r — 8.5 kpc) in Q. Three 
different halo models are shown with different values of 
the outer slope of 7 = 3, 4, 5. The inner slope is set in all 
cases to a = 1. The implied local density of dark matter 
is similar for each value of 7, and a larger scale radius is 
attained for steeper outer density slopes. This scaling of 
Tq with 7 is ultimately a reflection of the fact that the 
ratio of the escape velocity to the LSR velocity is bound 
to be 2.0 < w CS c/fisr ~ 3.5. 

Fig. [2] shows velocity distributions, vf(v), as a func- 
tion of 1 — v/v csc resulting from the numerical solution 
of the Eddington inversion formula for the three density 



models with outer slopes of 7 = 3, 4, 5 using the regions 
delineated in Fig. [TJ These fits are of dark matter halos 
with inner slopes of a = 1. The disk and bulge potential 
are included in the gravitational potential. The numer- 
ical evaluation of the Eddington formula extends over a 
wide enough range to ensure that the velocity distribu- 
tion function is convergent; in practice, this requires that 
the potential goes to zero at r > 100r s . Tangent curves 
show the analytic behavior of the outer slopes. The in- 
clusion of the disk component does not alter the behavior 
at small 1 — v/v csc . The bottom panel of Fig. [2] shows 
the error of using the ansatz velocity distribution func- 
tion of ([T]) or the King model instead of the full numerical 
solution of the Eddington equation, /e(v)- The error is 
defined as 



Error 



fk(v)-Mv) 



h(v) 



(22) 



The agreement between the power-law slope and the nu- 
merical calculation is the best for the density models with 
the steepest outer slopes, while numerical calculation for 
the NFW model with outer slope behavior of r~ 3 has a 
slightly steeper slope than k — 1.5. For the NFW, this 
is a reflection of the fact that the tail of the distribution 
is modified, so that as £ — ¥ 0, the distribution function 
scales as 



/(£) 



£1 



(-ln£) 



(23) 



rather than as a pure power-law in £ [2]. The behav- 
ior near £ = is steeper than a k = § power-law and 
/cnfw — 2.0 typically fits better numerically. 



B. Comparison to Simulated Galactic Halos 

High resolution simulations of Galactic halos with dark 
matter only have revealed a non-Maxwellian structure to 
velocity distributions [51 [7J. Several trends are seen in 
these simulations. The low velocity tail is more popu- 
lated than the best fitting Maxwellian distribution, while 
the peak is more depressed. At moderately high veloci- 
ties beyond the peak, the distribution is more populated 
than the best fitting Maxwellian. Further, broad bumps 
that represent dynamical features are prevalent in the ve- 
locity modulus distribution, though they are not present 
in the distribution of the individual velocity components. 
Although features in the distribution due to substructure 
are present at some locations in the halo, their contribu- 
tion to the velocity distribution averaged over all radial 
shells is typically sub-dominant relative to the broader 
bumps. The highest resolution simulation in the litera- 
ture that includes baryonic physics [8] does not show such 
features, and is better described by a Maxwellian distri- 
bution than the dark matter only simulations, though it 
still deviates at the high velocity tail. It remains to be 
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FIG. 3: Best-fit Maxwellian (dot-dashed), King (dashed), and Tsallis (dotted) models to the Via Lactea, Aquarius, GHALO, 
and Ling et al. [8] simulations (solid black). The first three are simulated with dark matter only; [8] simulates a galaxy with 
baryons. The green band shows the velocity distribution obtained by fitting ([T]) to the data for 1.5 < k < 3.5. The purple 
shaded regions enclose 68% of all measured distributions in each of the Aquarius, Via Lactea, and GHALO panels. 



seen whether this result is robust to changes in uncer- 
tain baryonic physics and scatter in formation histories 
between galaxies. 

The features of these velocity distributions clearly re- 
flect the formation process of each individual halo, and 
deviate from sphericity and isotropy in both the den- 
sity distributions and velocity ellipsoids. Of course, they 
also reflect the distribution of the individual dark matter 
halo that was simulated, and it is not yet known how this 
quantity in individual halos relates to the distribution of 
the velocity distribution for a larger sample of similar 
mass halos. Further, it is also assumed that the veloc- 
ity distributions from numerical simulations provide an 
unbiased tracer of the velocity distribution of the actual 
Milky Way dark halo. This latter point can be addressed 
by choosing simulated halos to have a mass and isola- 
tion criteria that resemble the Milky Way regime. These 
selection criteria could be further refined by demanding 
additional constraints, e.g., requiring a similar number 
of massive, Magellanic cloud-like satellites to the Milky 
Way [37] • In spite of these caveats, comparisons to N- 
body simulations are still critical for determining how 
well analytic arguments based upon the assumption of 



equilibrium match typical galaxies. 

Fig. [3] shows the best fits of the velocity distribution 
function from to the results of the Aquarius, Via 
Lactea II, and GHALO dark matter only simulations, 
and the simulation with baryons of Ling et al. J8J. The 
high energy tails are shown in the main figures in each of 
the four panels while the insets show the fits to the full 
velocity distribution function. The shaded bands in this 
figure represent slopes in the range k — [1.5,3.5], corre- 
sponding to outer slopes between 7 = [3.0, 5.0]. For each 
value of k within this band, the parameter vq has been 
fit to the distribution, and v csc is set to the value de- 
termined by each of the respective simulations. Implicit 
in this assumption is that the highest velocity particle 
is a faithful tracer of ^ eS c, rather than a contaminating 
particle that is not bound to the halo. 

The best-fit SHM, King, and Tsallis distributions are 
shown for comparison in all panels of Fig. [3] In all pan- 
els, the SHM does not go to zero smoothly at the escape 
velocity, indicating that this distribution function is un- 
physical in the tail. By construction, the velocity distri- 
bution in ([!]) smoothly joins to a Maxwellian distribution 
for any value of k, with the differences manifest only near 



r 



the tails. 

The fits shown in Fig. [3] indicate that the velocity dis- 
tribution in ([I]) captures the high velocity tail of the sim- 
ulated distribution well, especially in comparison to the 
King, SHM, and Tsallis distributions. This is in spite of 
the fact that the simulations are clearly anisotropic and 
non-spherical, and that the outer regimes of the halo may 
not have fully equilibrated. The agreement demonstrated 
in Fig. [3] is promising, though larger samples of dark 
matter halos at similar resolution are needed to verify 
the high velocity suppression that ([I]) predicts. Consis- 
tent with previous findings, the dark matter only simula- 
tions are more 'flat-topped' than all of the models, while 
there are more particles in the low regime of the veloc- 
ity distribution. The enhancement of low velocity parti- 
cles may result from radially-biassed anisotropy, as dis- 
cussed below. The Ling et al. [8] distribution best resem- 
bles a Maxwellian profile at low velocities. This agree- 
ment might be attributed to either the effect of baryonic 
physics, which has been shown to result in more isotropic 
inner halos [38 , or to the lower resolution of dark matter 
simulations that include baryons. 

A final point regarding the comparison to the simu- 
lations is that, since the tail appears well described by 
k = [1.5,3.5], it is plausible to conjecture that baryons 
do not appear to greatly affect the shape of the velocity 
distribution near the escape velocity. One way to un- 
derstand this is that since highly energetic particles with 
large peri-center orbits shape the tail of the distributions, 
these dark matter particles will be outside the halo's scale 
radius for a large fraction of their orbits. The scale ra- 
dius ranges between r s ~ [10, 60] kpc (see Fig. [I]) for the 
models considered, and the Galactic disk is located in 
the regime where the density scales as the inner slope, 
a. The inner slope should not significantly perturb the 
orbits of the halo's most energetic particles that make up 
the high velocity tail. Thus, the regime most likely to be 
affected by the formation of the disk are the orbits in the 
the inner core of the halo, where the power-law depends 



on a. 



IV. DIRECT DETECTION 

The velocity distribution function in ([I]) reproduces 
the results from N-body simulations and exhibits non- 
Gaussian behavior near the escape velocity. This section 
shows that the suppression in the high energy tail has 
important implications for the interpretation of direct 
detection experiments. 

The scattering rate of a dark matter particle off a tar- 
get nucleus depends on the velocity-averaged differential 
cross section and is thus sensitive to the distribution func- 
tion. The detection rate per unit detector mass is 



dR k 
dEn 



Pa 



d 3 v f k {v + v E {t))v 



da 
dEn 



(24) 
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FIG. 4: Fractional change in the differential dark matter de- 
tection rate as a function of the minimum velocity for k = 1.5 
(red), 2.5 (yellow), and 3.5 (green), compared to a King veloc- 
ity distribution function with k = 1. In this figure, vo = 220 
km/s, Vcsc = 550 km/s, and the Earth's velocity is taken at 
~ June 2. 



dark matter mass, hn is the reduced mass of the nucleus- 
dark matter system, Er is the nuclear recoil energy, po 
is the local dark matter density, and fk{v + vsit)) is the 
velocity distribution function in the Earth's rest frame 
with a power-law index, k |211 I39j . The spin- independent 
differential cross section is parameterized as 



da m^a n 

dER = 



(f p Z + f n (A-Z)) 2 \F N (q 2 )\ 2 , (25) 



where a n is the dark matter-nucleon cross section at zero 
momentum transfer, fj, n is the dark matter-nucleon re- 
duced mass, and q 2 = Im^E^ is the momentum trans- 
fer- fp.n parameterize the coupling to the proton and 
neutron and are typically set to unity. FN(q 2 ) is the nu- 
clear form factor, for which the Helm/Lewin-Smith form 
is used [2"T] . 

The relevant quantity when comparing direct detection 
rates for velocity distributions with different spectral in- 
dices is 



9k(v B 



dv v f k (v + v E (t)). 



(26) 



where m^r is the mass of the target nucleus, rridm is the 



The ratio gk/91 illustrates how the new ansatz compares 
with the more commonly used King model. In the limit 
where v m - m — > 0, the ratio gk/91 approaches unity for 
properly normalized distributions. However, this ratio 
deviates from unity when the minimum scattering ve- 
locity is ~ 60% of the escape velocity, as illustrated in 
Fig. |U As the minimum velocity approaches the max- 
imum dark matter velocity, v csc + Ve(£), the scattering 
rate is strongly suppressed. This suppression can be as 
large as two orders of magnitude for an outer slope of 
7 = 4, or five orders of magnitude for 7 = 5. 

Dark matter with particularly large minimum scatter- 
ing velocities will be most sensitive to increases in the 
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FIG. 5: Fractional change in the scattering rate for k = 2.5 compared to As = 1 for an elastically scattering dark matter with 
vq — 220 and v csc = 550 km/s. The plot on the left illustrates the dependence of the scattering rate on the dark matter mass for 
DAMA-Na (dot-dashed) g2], CoGeNT (solid gray) [23], XENON (dashed) g3H5], CDMS (solid black) gBJ, and CRESST-W 
(dotted) [47] for threshold energies of 7.5, 1.7, 5, 10, and 10 keV, respectively. The plot on the right illustrates the dependence 
of the scattering rate on the threshold energy for an 8 GeV dark matter scattering elastically off a Xe (dashed), Ge (solid), and 
Na (dot-dashed) target. In both plots, the Earth's velocity was taken at ~ June 2. 



spectral index of the velocity distribution function. The 
minimum velocity depends on the kinematics of the scat- 
tering event and is given by 

v 1 f rn N E R | a 

y/2m N E R \ ' 

where 5 is the mass splitting between the initial and fi- 
nal dark matter states if dark matter scattering is dom- 
inantly inelastic; elastic scattering corresponds to <5 = 0. 
Figure [4] shows the effect of the power-law index k on the 
direct detection rate, as a function of v m i n . 

Inelastic dark matter [40] is sensitive to the spectral 
index of the distribution function because it has a larger 
Vmin- The minimum scattering velocity differs between 
experiments with different energy thresholds or target 
nuclei. As a result, a change in the spectral index of 
the velocity distribution affects the limits for each ex- 
periment differently. For example, an experiment with 
a heavier target nucleus or a lower threshold energy will 
have a fairly large value for w m i n and a greater sensitivity 
to the velocity distribution function. 

Light elastic dark matter is also sensitive to the high 
velocity tail of Galactic dark matter gT] . Fig. [5] illus- 
trates the suppression of the scattering rate Rk (inte- 
grated over recoil energy) relative to that for the King 
model. For example, if md m = 8 GeV, a Xe target with 
5 keV threshold energy has R2.5/R1 — 0.18. This ratio 
is more suppressed as one goes to heavier target nuclei 
and/or larger threshold energy. For instance, a Ge tar- 
get with E tn = 10keV is suppressed relative a Ge target 
with a E tn — 2 keV by a factor of 25. In the case of light 
target nuclei (i.e., Na), or low threshold energies (i.e., 
< 2 keV) , the effect of increasing the spectral index is 
much less significant, with R2.5/R1 ~ 1. The relative 
changes in the rates have profound consequences when 



comparing a potential signal in one experiment with the 
result of another experiment with a different target nu- 
cleus or threshold energy. 

V. DISCUSSION AND CONCLUSION 

This paper has considered a cosmologically-motivated 
class of double power-law dark matter density profiles 
and has examined the resulting velocity distribution, as- 
suming isotropy and spherical symmetry. This article 
has specifically focused on the behavior of the velocity 
distribution at low binding energies, near the escape ve- 
locity. The primary result is that the velocity distribu- 
tion function is well-modeled by the function presented 
in ([I]), where the power-law index ,fc, is in the range 
~ 2. This is a steeper fall-off than what is expected 
for a Maxwell-Boltzmann distribution, which has served 
as the canonical model for dark matter velocity distribu- 
tions in the recent literature. The new ansatz in ([lj af- 
fects the interpretation of limits set by direct dark matter 
detection experiments. Limits on low mass dark matter, 
mdm < 40 GeV, are most strongly affected, especially for 
experiments with heavy target nuclei or large threshold 
energies. 

The analysis presented in this article does not address 
the effects of anisotropy and deviations from spherical 
symmetry, which are expected to be more serious in the 
outskirts of the profile most relevant for the high velocity 
tail [IS] . For instance, anisotropy could be responsible for 
a <C 2. For anisotropic distributions, there is no guar- 
antee that a given density distribution will correspond 
to a unique distribution function. For example, if the 
distribution function is deformed to L~ 2/3 f(£), the dis- 
tribution function depends only on a single derivative of 
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the density with respect to the potential for radial orbits 
of y9 = | [H]. The corresponding velocity distribution 
scales as f(v) oc v at low velocities, which is broadly con- 
sistent with the behavior in the dark matter only simula- 
tions of Fig. [3] In comparison, the low velocity behavior 
of the Galactic halo for a circularly biased system with 
/3 = — \ is f(v) oc v 3 . Thus, the low energy tail may 
indicate the presence of radial orbits in the halos, which 
would be consistent with the velocity anisotropy distri- 
bution in numerical simulations [49] . 

The interplay between density profiles of dark matter 
halos and their corresponding velocity distribution func- 
tions was developed in this paper, but requires further 
study. Exploring equilibrium dark matter phase space 
distributions, particularly anisotropic distributions, may 
give further insight into the behavior of N-body simu- 
lations and may explain the inner slope of p(r). These 
effects are unlikely to be directly observable in early di- 
rect detection experiments, but may be seen in future 
directional dark matter detection experiments |50| . 
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